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fvl Abstract 

J>~> The discovery of the " Standard Model Scalar" (Brout-Englert-Higgs boson) opens a new 

^2 field of research, namely the structure of the scalar sector. Numerous extensions exist, and 

^^ imply extra particles, including additional scalars. Two main alleys are open to investigate 

^_ such deviations from the minimal Standard Model: precision measurements could indicate a 

fvj deviation from the usual expectations, or a direct discovery of new scalar partners (or other 

particles) would establish an alternative. In this short note, we concentrate on a very simple 
model, for which the respective reaches of the two approaches can be compared. This note 
also provides a strong incentive to pursue searches for extra "Standard Model Scalar- like" 
particles in the whole available energy spectrum. 
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CN 1 Introduction 

> 

A new scalar has been found at LHC, with a mass around 125 GeV. To check that this is indeed 
^r^ the minimal version of the Brout-Englert-Higgs [I] [2] boson in the Standard Model [3], we can 

•O either rely on increasingly precise measurements comparing branching ratios and production 

■^ rates, or satisfy ourselves that no further scalar structure is found. Both are, of course, open- 

ended tasks, as extensions of the Standard Model scalar structure can in principle be arbitrarily 
close to the basic version. 

To compare the approaches, we rely on a very simple model. Its initial version dates back 
more than 30 years f5], but it is found as an ingredient in a number of more elaborate construc- 
^ tions, including the Next to Minimal Supersymmetric Standard Model (NMSSM) [6j. 

C^ The principle of the model is very simple, and we will describe it quickly, before providing 

the corresponding Lagrangian: add a (real) singlet scalar S boson to the minimal Standard 
Model doublet H. As a singlet under the gauge group, it does not interact with the gauge 
bosons, and has no direct coupling to the fermions. Its only couplings to known particles are 
through the scalar. After symmetry breaking, the usual "Standard Model Scalar" (BEH boson) 
IL mixes with the neutral singlet. We end up with two neutral scalars: 

Hi = Hcosa + Ssina (1) 

H2 = —Hsm.a + Scosa. (2) 



Since the singlet S possesses no interaction of its own, the Hi and H2 states interact only 
by their H component. We thus end up, at this stage, with two mass states (2 "peaks"), rrii 
and m,2, and will by convention take rrii = 125 GeV, while m2 can be heavier or lighter. 

Note that in the balance between precision measurements and direct discovery, it should 
be kept in mind that the present example is a kind of "worst case" for the latter approach. 
Indeed, the particle added to the Standard Model exhibits no new interaction of its own, and 
only interacts through mixing. More general extensions would involve new interactions, which 
could boost production of the new components. 

For simplicity, we will use SMS for the canonical "Standard Model Scalar". In terms of 
production, the Hi and H2 particles are produced each exactly like a SMS of the corresponding 
mass, but with factors cos^ a and sin a respectively. For their decay, the branching ratios are 
identical to those of a hypothetical SMS of the same mass mi or m2- Of course, the total width 
of the state, like the production rate, are weighted by cos^ a and sin^ a respectively, but this 
does not affect directly the observation rates in the various channels (once produced, unless a is 
vanishingly small, both particles decay quickly). One indirect effect might appear in the width 
of the peak, when this is observed - this would facilitate the detection of heavy scalars, but 
probably requires the adaptation of the search. 

It is thus fairly straightforward to exploit the existing searches for a SMS at various energies 
(and the corresponding bounds) to explore this model. For the time being, we will consider the 
parameters mi = 125 GeV, m.2 and a as independent (see more details below). 

We can thus present the existing data (possibly channel by channel, but preferably com- 
binations, to achieve the best sensitivity in a large energy range) in different ways within the 
context of the model. The simplest graph would provide the limits (or discovery) for the whole 
explored mass range, and the ratio (^x_b_r) ''^ ^^ ordinates represents the limits on cos^ a and 
sin^ a respectively. 

In order to compare the "reach" of the precision measurement techniques however, we prefer 
to plot the limits obtained on cos^ a and sin'^ a as ordinates and abscissae respectively. The 
minimal model considered here is then represented by the line cos^ a+sin^ a = 1, and it becomes 
easy to check which constraint (departure from expected rate at 125 GeV or limit on an extra 
peak at another energy) is the most restricting. Of course, upper limits on sin^ a will depend 
on the energy range considered, and a series of exclusion lines will need to be drawn. 

We have performed the exercise in Figfll using currently available (public) CMS data [7]. 
Unfortunately the confidence levels of published direct searches (95% C.L.) and production rates 
(1(t) are different. We have plotted them as such (since we don't know the exact x^), but one 
should keep in mind that the horizontal band corresponding to the production rate should be 
considerably broadened for a fair comparison to the direct searches. The result is nevertheless 
quite interesting, as it shows that, except for the lowest (below 125 GeV) and highest (above 
650 GeV) values of the masses, the direct search is most constraining, despite being disfavoured 
by the specific model. We hope that this comparison can serve to plot the expectations of future 
LHC runs, upgrades and alternative machines. 

An alternate way to plot the constraints of this simplistic model consists in superimposing 
the upper bound on sin^ a obtained from the 125 GeV scalar production rate as an horizontal 
line on the usual SM-like search plot z'^xBR) ^® mass of the extra scalar; it provides a more 
synoptic view, but insists less on the complementarity of the approaches in narrowing the model 
parameters. 
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Figure 1: The proposed benchmark compared to current pubhc CMS data (HCP 2012) [7j. 
The extra singlet benchmark is constrained to be on the blue oblique line sin^ a + cos^ a = 1. 
The ordinates correspond to the precision measurement of the production of the 125 GeV SMS 
(ax BR) ' ^^^ pink hatched area represents the current la confidence interval. The abscissa 
gives current production limits (this time at the more stringent 95% C.L.) for SM-like scalars 
of different masses (110 GeV, 160 GeV, 200 GeV, 300 GeV, 600 GeV and 650 GeV). 



2 Simple models 

We have treated this far the two main parameters (the mixing angle a and the second scalar 
mass 7712) as independent. We will show below that some dependence exists (as implied by 
decoupling requirements), but that it can be safely ignored at this stage. For this purpose, we 
now present an explicit form of the (scalar part of the) LagrangiairJ - to be further referred to 
as Extra Singlet Model (ESM). In [5], the authors introduce the following Lagrangian (<^ and x 
are the doublet and the singlet respectively )rl 
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(3) 



It is an easy task to check that |$p = fi/2 and x = /f/2/2 minimize the potential. Then we 
can identify /i with v ~ 250 GeV, the usual VEV of the SMS field, linked to the W mass. The 
mass eigenvalues are then: 
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^A detailed study of the scalar potential stability in a related model can be found in [4]. 
^The present normalization differs by an (arbitrary) factor of 2 from the one of reference [5]. 



(4) 



with A3 = Ai + A2- One can then show that the mixing angle is given by (both for mi > 7112 
and m,2 > mi): 

A3 - [mi/vf 



sin Q 



(5) 



{m2/v) — {mi/v) 

As expected, decouphng is achieved in the hmit where m2 becomes much larger than mi with 
bounded A3, as the mixing vanishes accordingly. 
The alternative notation: 
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(6) 



shows that, once m-i is fixed, each mixing a is associated with a straight line in the plane 
{{m,2/v)'^, A3). Figl2 shows the two extreme lines (for mi/v = 1/2): sin^ q = is the horizontal 
red line, while sin^ = 1 is the oblique red line. The whole set of intermediate values for 
the mixing corresponds to lines that lie in between (red hatched region). They all cross when 
m2 = nil. It can be shown that the two extreme lines must be excluded from the set of physical 
parameters. One can understand it easily: since in this simple model the mixing is proportional 
to one of the diagonal elements of the mass matrix, then, if there is no mixing (extreme lines 
cases), one cannot find a solution for a non zero m,2- 
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Figure 2: Allowed region of parameters for the model (JSJ). The red hatched region gives the 
physical values (the border are not physical, but are included with the additional A4 coupling (see 
text for details)). Even for modest A3 most of the (a,rn-2) parameter space remains accessible. 

The only dimensionless parameters that are relevant in the discussion are Ai and A2. Stability 
forces them to be positive, and then the perturbative regime is insured while A3 remains small. 
The meaning of small is quite subjective here (as it is difficult to judge the convergence of the 



series), so we simply provide a plot (FigM linking A3, a and 1712- Even for very conservative 
values of A3, a large range of parameters (0,7712) is open. 

Just for completeness, the extreme values of sin^ a = and sin^ a = 1, excluded in the basic 
Lagrangian, can be recovered by adding an extra term: 




^V = -^ r-Urr , (7) 



with a new dimensionless parameter A4. This term is chosen such that the minimum of the 
potential is unchanged. 

In the large 7712 limit, decoupling is achieved as expected, but in a relatively slow manner 
(e.g. with A3 < 2, a mixing of 30% it still allowed for 7712 ^ 1 TeV). Indeed a much faster 
decoupling is usually found in comparison, for instance in the Two Doublets Model (2DM). The 
most general potential for such a model can be parametrized as: 



i2 rf>t,f>. I ^2 ^t 



t2 ^t. 



V = mli(^[(^i + m22^2^2- 7ri^2^i*^2 + h.c. 

+ y ($I$l)2 + ^{^\^2? + A3($I^l)(^$2) + \a{^\^2){^\^i) (8) 

+ |^($1$2)' + [A6(^I$i) + h{^\^2)\ $1^2 + h.C.j . 

In general 777|2, A5, Ae and A7 can be complex, with possible CP-violating effects. We will ignore 
this possibility here, as it is unrelated to our problem. If the mass matrix possesses at least one 
negative eigenvalue, the scalar fields will develop a VEV. Imposing CP invariance and U{1)em 
gauge symmetry, it can be written (<^j) = (0 Vi/^/2)'^ , with Uj G M and uf + w^ = f^, linked to 
the W mass. Then, it is always possible to choose a basis wher^ui = and V2 = v. 

It is then easy to derive the following relations for mi =125 GeV, 7712 and a, the masses and 
mixing of the two neutral CP-even scalars (and also ruA and mH± the masses of the neutral 
CP-odd and charged scalars that are not eaten by gauge bosons): 
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where \ij,„k = \i + Xj + ... + A^. 





m\ + Asf^ = 777^± + -A45U^ j (9) 



(10) 

(11) 



As in the ESM, there are two dimensionfull parameters v and 77711 that come into the game. 
While the first one is fixed by the W mass, the other one can be used to push all the masses 
(expect 7771 ) to a huge value (they all become almost equal, because the Aj are bounded by 
perturbative arguments). The relation (10) is the equivalent of ([6]) in the ESM, and gives 



essentially the same constraint on the mixing: when 7772 becomes too large, a small sin a must 
compensate. 



The difference in the 2DM is the presence of the new relation (11). In the ESM, the off- 
diagonal term of the mass matrix contained a product of the two dimensionfull paramaters v 
and /2 rather than only v"^ like here. Then, when /2 is increased to raise 7772, this term "follows" 



^ When couplings to matter are introduced, the choice of the basis is no longer arbitrary. However, this choice 
permits to draw more easily a parallel with the previous model and we expect our conclusions to remain valid in 
any basis. This will be supported by a comparison with a more general result from [8]. 



and we find no new constraint. On the contrary in the 2DM the condition (11) introduces a 
new constraint at higher m2 as shown in FigJ3J 

In the 2DM model, this result was found in the same context, but in the more general case 
where ui 7^ [8]. The author showed that in the decoupling limit (7^2 ~ tua ~ 'm'H± ^ mi): 
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(12) 



where tan/3 = V2lv\, and m,^., w-t and rriD are some functions of the CP-even mass matrix 

(13) 



elements (see eqns (3.5) in ^). In our limit /3 — )■ 7r/2 and (12) turns into: 

2 



sin a 



A7 
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We see on Fig J3] that this line gives the same constraint as (11), in the regime 7712 S> ra\ 




m^ [TeV] 

Figure 3: Constraints on the mixing sin a as a function of m,2 in the present ESM case (upper 
curve, in red) and in the 2DM models (upper limit in thin black, approximation ( 13 ) in green- 



dashed) for A2,A7 < 10. 



3 Conclusions 

Extensions of the Standard Model (notably through extra scalars) can be studied either by 
searching for a deviation in the production rate of the 125 GeV Standard Model Scalar (Brout- 
Englert-Higgs scalar), or by direct search for additional states. Taking the Extra Singlet Model 
as an example, we see that direct detection is currently much more sensitive in most of the 
accessible mass range (125-600 GeV). 

Taking perturbative unitarity as a guide, a full coverage of the parameter space allowed in 
the future by a few percent departure in production rate would require to extend the detection 
sensitivity to the TeV range. The model considered here offers simplicity as a benchmark (the 
extra state has the properties expected from a " Standard Model Scalar" of the same mass, so 



6 



existing analysis can be directly used) . The fact that the extra component is a singlet (without 
its own direct interactions to known matter) makes direct detection more difficult, while the 
maximal mixing allowed by perturbative unitarity is rather large (at least compared to the 2 
Doublets Model). 
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